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Abstract - The most salient features of the surface structure and bulk behavior of Ni-Pd alloys 
have been studied using the BFS method for alloys. Large-scale atomistic simulations were per- 
formed to investigate surface segregation profiles as a function of temperature, crystal face, and 
composition. Pd enrichment of the first layer was observed in (111) and (100) surfaces, and 
enrichment of the top two layers occurred for (110) surfaces. In all cases, the segregation profile 
shows alternate planes enriched and depleted in Pd. In addition, the phase structure of bulk Ni-Pd 
alloys as a function of temperature and composition was studied. A weak ordering tendency was 
observed at low temperatures, which helps explain the compositional oscillations in the segrega- 
tion profiles. Finally, based on atom-by-atom static energy calculations, a comprehensive explana- 
tion for the observed surface and bulk features will be presented in terms of competing chemical 
and strain energy effects. 
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1. INTRODUCTION 


A complete understanding of the surfaces of bimetallic systems offers significant challenges 
because the structure, composition, and thus properties are usually quite different than that 
expected from the bulk [1], Because surfaces have an extra degree of freedom compared to the 
bulk, a number of phenomena unique to surfaces may be observed, such as surface segregation, 
surface ordering, and reconstruction. Better knowledge, understanding, and eventual exploitation 
of these phenomena can lead to the development of mixed-metal compounds with potential wide- 
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spread application in the areas of electrochemistry, microelectronics, and especially catalysis, 
where mixed metal systems have been found to be much more effective than pure metal surfaces 
[1,2]. The Ni-Pd system is one such example, where surfaces of the binary alloy are more efficient 
than the surfaces of either of the pure metals at hydrogenation reactions [3-5] and therefore have 
been studied extensively both experimentally [3-9] and theoretically [10-12], In spite of the 
amount of work performed on the Ni-Pd system, there are still questions concerning the surface 
and near-surface structure in this bimetallic system and the existence of ordered phases. 

Ni(Pd) alloys are an example of a dilute alloy system for which very high segregation of the 
noble metal to the surface occurs. In fact, it takes very little Pd in solution to significandy enrich 
the surface, with polycrystalline Ni alloys containing 1 and 5 at.% Pd consisting of surfaces 
enriched with 20 and 50% Pd, respectively [6]. In the study by Miegge et al. [6], the overall 
enrichment did not seem to extend beyond the second atomic layer and approached the bulk com- 
position by the third layer. Additional analysis of this same Ni-5Pd alloy confirmed that only the 
surface layer itself was enriched [3], In fact, significant Pd-enrichment of essentially just the sur- 
face layer is consistent with all polycrystalline alloy studies [6,7,9]. The driving force for this 
behavior can be rationalized by the fact that Pd has a lower surface tension than Ni, which in addi- 
tion to its larger atomic radius, results in expulsion of the noble metal from the bulk and its segre- 
gation to the surface [13]. 

Studies of single crystal surfaces are consistent with the polycrystalline data in that they also 
show significant surface enrichment in Pd [4,8]. These single crystal studies also revealed addi- 
tional structural detail, including oscillatory composition profiles in the near-surface region. The 
first to observe this phenomenon was Derry et al. [8] in a single crystal (100) Ni-50Pd sample. As 
expected, the surface layer was enriched in Pd (80%). But the first layer below the surface was 
observed to be 100%Ni, and the second layer below was enriched once again in Pd (64%), though 
less so than the surface layer. 

The observed oscillatory behavior is not easily explained in terms of simple segregation mod- 
els, especially since Ni-Pd is considered a prototypical solid solution alloy over the entire binary 
alloy composition range with no reported ordered intermetallic phases [14]. Consequently, obser- 
vation of an oscillatory segregation profile in the dynamical LEED study of (100) Ni-50Pd by 
Derry et al. [8] is quite unusual [15] and until now has not been fully explained. Furthermore, 
recent analysis has shown that in addition to the results of Derry et al. [8] where the (100) surface 
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layer is rich in Pd and the second layer is essentially pure Ni, results by Michel et al. [4] on (1 10) 
surfaces strongly indicate that the first two layers are highly enriched in Pd. No single explanation 
has yet been presented to explain whether this difference in behavior as a function of surface ori- 
entation is consistent and to what extent the known surface behavior is compatible with the behav- 
ior of bulk Ni-Pd alloys. 

In light of these questions, we have used the Bozzolo-Ferrante-Smith (BFS) method for alloys 
[16] to study Ni-Pd alloy surfaces, including layer-by-layer compositional analysis of segregation 
behavior for different crystal faces. Similar analyses were performed for determination of the 
phase structure of bulk Ni-Pd alloys. From this analytical study we can provide a consistent expla- 
nation for the previous experimental findings [4,8] in terms of surface energy arguments, strain 
effects, and equally important in this case, bulk alloy behavior. 

2. THE BFS METHOD 

The basic features (both surface and bulk) of the Ni-Pd system can be described within the 
framework of a simple modeling approach grounded in the BFS method for alloys. The BFS 
method [16] is based on the concept that the energy of formation of a given atomic configuration 
is the superposition of the individual atomic contributions, A H = £ e, 7 ! Each contribution e, r is the 
sum of two terms: a strain energy, e,^, computed in the actual lattice as if every neighbor of the 
atom i was of the same atomic species i, and a chemical energy, e, , computed as if every neigh- 
bor of the atom i was in an equilibrium lattice site of a crystal of species i, but retaining its actual 
chemical identity. The computation of £,' < \ using equivalent crystal theory (ECT) [17], involves 
three pure elemental properties for atoms of species i: cohesive energy, lattice parameter, and bulk 
modulus. The chemical energy, £j C , includes two BFS perturbative parameters, A NiPd = -0.0396 
A ' 1 and A PdNi = -0.0478 A' 1 , computed by fitting the BFS calculations to the experimental values 
of the heat of solution in the dilute limit [18]. A reference chemical energy, e, Co , is also included 
to insure a complete decoupling of structural and chemical features. Finally, the strain and chemi- 
cal energies are linked with a coupling function g, which ensures the correct volume dependence 
of the BFS chemical energy contribution. Therefore, the contribution of atom i to the energy of 
formation of the system is given by 
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Large scale simulations were performed using the Monte Carlo - Nearest Neighbor (MC-NN) 
approach described in Ref. 19, which is a modified version of the traditional Monte Carlo - 
Metropolis algorithm. An initial, random, structure is generated by assigning Ni and Pd atoms to 
rigid (i.e., no individual relaxations) lattice sites in the computational cell with periodic boundary 
conditions in two directions for surfaces (three directions for bulk cells). Segregation to surfaces 
or bulk ordering are then simulated by the exchange of randomly chosen pairs of atoms of differ- 
ent atomic species throughout the cell which only occupy nearest-neighbor (NN) sites. The differ- 
ence in energy of formation of the cell before and after the exchange, A E, is used to determine the 
likelihood that this exchange can take place. BFS was used to calculate the energy at each step. If 
A E < 0, the exchange is readily accepted. If not, it is assigned a probability exp(A E/Ejj), where 
the available thermal energy, Ej 7 ,, is given by 



( 2 ) 


where 0 D is the Debye temperature [20] and k is Boltzman’s constant. As such, the MC-NN 
approach does not rely on a rigorous statistical mechanics foundation as traditional Monte Carlo- 
Metropolis does. It is, at best, a simple but approximate way to model the evolution of the system 
[19]. 

Table 1 lists the necessary parameters for applying the BFS method to the Ni-Pd system. We 
refer the reader to Ref. 16 for a detailed discussion of the BFS method, its definitions, operational 
equations and their implementation. All calculations in this work were performed using the adw- 
Tools software package [21]. 

While the methodology allows for a detailed analysis of individual and collective atomic 
relaxations, no such effects have been included in this initial study. The purpose of the paper is to 
investigate the balance between strain and chemical effects as driving forces for segregation and 
ordering, thus focusing on a simple but unified description of bulk and surface behavior. 
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3. RESULTS AND DISCUSSION 


3. 1 Surface and near-surface behavior 

The computed segregation profiles, determined from MC-NN simulations for Ni-5Pd and Ni- 
50Pd alloys, as a function of low-index crystal face, are shown in Fig.l. Each figure includes 
results for T=300 and 1000 K. Two features are common to all the cases shown: first, the segrega- 
tion of Pd to the surface layer (top two layers in the case of a (110) face), and second, damped 
oscillations in the concentration profile for the following subsurface layers. These oscillations, 
which are usually associated with ordering tendencies in the bulk alloy [15], are more noticeable 
for the T=300 K profiles. These two basic features are independent of concentration, as all the 
profiles share a similar Pd-surface enrichment and oscillatory segregation profile near the surface. 
The segregation profiles also exhibit rapidly decreasing compositional oscillations in the subse- 
quent subsurface layers, approaching the bulk composition after several layers, depending on the 
bulk concentration of Pd. 

While these simulations correctly determine not only the gross behavior but also fine details of 
segregation in the Ni-Pd system as described in Refs. 4 and 8. Further details on the correspond- 
ing driving mechanisms for these observations can also be obtained within the framework of a 
simple modeling approach. An atom-by-atom energy analysis of a single Pd atom on or near a Ni 
surface can be enough to explain the observed differences between the segregation trends in 
Ni(110) and other surfaces. Having shown in Fig. 1 that the (100) and (111) surfaces display the 
same segregation trends, but different from that of the (110) surface, we limit the comparison to a 
(100) and a (1 10) surface. Fig. 2 displays a side view of a ( hkl ) slab, with the Pd atom located in 
different planes. We will denote with Pd(L//ifc/) the configuration where the Pd atom is located in 
layer L (L = S, lb, 2b, representing the surface, one atomic layer below the surface, and two lay- 
ers below the surface, respectively) parallel to a (hkl) surface. The role of each ‘active’ atom (i.e., 
the Pd atom and all the Ni atoms noticeably affected by the presence of the Pd atom) is quantified 
by its individual contribution (e,- 7 ) to the energy of formation (AH) of the cell and is shown in the 
corresponding table, along with the BFS strain and chemical contributions that make up e? (Eq. 
1). The net difference in energy arising from each active atom with respect to an equivalent atom 
in a pure Ni surface, AE/atom, and the number of identical atoms ( m ) with the same energy in the 
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cell are also shown. A reduction in energy occurs when AE/atom is negative, since it represents a 
decrease in energy with respect to the contribution of that atom to AH in the VdQL/hkl) cell com- 
pared to the contribution of a Ni atom in the same site in a perfect Ni(100) or Ni(l 10) surface slab 
(i.e., no Pd present). Whether the entire system is more stable or not depends on the total contribu- 
tion to the energy by all affected atoms, which is determined by adding the values in the AE/atom 
column (once these values are multiplied by the multiplicity factor m). 

For example, Fig. 2.a shows a Pd atom in the (100) surface plane (Pd(S/100)). The Pd atom 
contributes 0.4597 eV to AH, which amounts to a decrease of 0.7286 eV compared to an equiva- 
lent site in a pure Ni(100) surface. Each of the four surface Ni nearest-neighbors (NN) of the Pd 
atom contribute 1.1693 eY to AH, resulting in a total decrease of 0.0760 eV with respect to a pure 
Ni(100) surface slab. Other surface Ni atoms, next-nearest-neighbors (NNN) of Pd or beyond, 
introduce negligible changes in energy, -0.0001 eV/atom or less, compared to sites in a Ni(100) 
surface. Finally, each of the four Ni NN in the lb layer contribute -0.0027 eV/atom to AH, intro- 
ducing a total decrease in energy of 0.0108 eV. Ni atoms in the 2b layer are barely affected by the 
presence of Pd(S), thus having a negligible contribution to AH. The analysis indicates that all the 
Pd NN Ni atoms provide favorable contributions to AH resulting in a lower formation energy, but 
the changes are an order of magnitude smaller than that due to the Pd atom itself. The presence of 
the one Pd surface atom lowers the energy of the cell by 0.7942 eV. This decrease in energy is a 
result both of the lower surface energy of Pd and the favorable bonding between Ni and Pd atoms 
(as denoted by the negative contribution to AH from the affected Ni atoms). 

Fig. 2.b shows the Pd atom in the first plane below the (100) surface (without relaxation) 
(Pd(lb/100)). The coordination of the Pd atom is nearly bulk-like (with the exception of one miss- 
ing NNN) and therefore the strain due to the larger Pd atom is too large to introduce a net decrease 
in energy to the system. The favorable bonding, however, does affect the surrounding Ni atoms 
introducing a slight decrease in their energy, but not large enough to offset the strain effect due to 
the embedded Pd atom, resulting in a net gain in energy for the cell of 0.1426 eV. A similar situa- 
tion arises for the Pd atom located in the second plane below the surface (Pd(2b/100), as shown in 
Fig. 2.c. 

A different picture emerges from a similar analysis of a Ni slab with a (1 10) termination. Fig. 
2.d shows the case of a Pd atom in a (1 10) surface site, (Pd(S/l 10)). As in Pd(S/100), most of the 
energy gain arises from the lower surface energy of Pd, thus guaranteeing Pd enrichment of the 
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top layer. The lower coordination of a Pd atom in a lb site (Fig. 2.e), however, leads to an interest- 
ing result: the Pd atom, surrounded by just 10 NN, finds itself in an environment of almost negli- 
gible BFS strain. In other words, the strain energy term, £ { s , as calculated with ECT [16,17], is 
essentially zero. This fact, added to the favorable nature of the Ni-Pd bonds that are created, low- 
ering the contributions to AH of the surrounding Ni atoms, leads to an overall decrease in energy 
(-0.3124 eV) for the cell, although not as large as the one computed for the Pd(S/l 10) case. It is, 
however, large enough to favor an enrichment of the lb layer, as observed experimentally [4]. It 
should be noted that while the insertion of Pd in the lb (1 10) surface layer results in energy gains 
due to a reduced strain, this fact alone is not sufficient to explain the energetically favored Pd(lb/ 
1 10) configuration which, again, arises from a combination of low strain and favorable Ni-Pd 
bonds. For Pd(2b/1 10), however, as in the case of Pd(lb/100) and Pd(2b/100), the strain is too 
large to be offset by a small chemical energy gain, thus resulting in segregation to the upper layers 
and depletion of Pd in 2b (110). 

This analysis, however, is not sufficient to explain the oscillatory nature of the segregation 
profiles immediately below the surface. But the favorable Ni-Pd bonds observed in the previous 
analysis does suggest a mechanism for the alternating layer eririchment/depletion of Pd in the 
near-surface layers. Consider the (100) surface where the top layer is enriched with Pd. Fig. 3 
demonstrates two (idealized) possibilities: a single Pd atom migrates to a site in the lb plane (Fig. 
3.a) or the 2b plane (Fig. 3.b). Following the notation used in Fig. 2, ‘active’ atoms are identified 
by their individual contributions to AH and by the net change in energy with respect to a reference 
state which, in the examples in Fig. 3, corresponds to a pure Ni(100) or (1 10) slab with a complete 
Pd surface layer (i.e., no Pd in layers below the surface) and an additional Pd atom somewhere 
deep in the bulk. Let A(L) denote an atom of species A (A = Ni, Pd) in layer L (L = S, lb, 2b). In 
the first case (Fig. 3. a), the Pd(lb) atom is in a high-strain site (1.0795 eV) due to the almost-full 
coordination of a lb site in a (100) slab, resulting in a large positive contribution to the energy of 
the system (0.7697 eV). The Ni atoms in the lb layer that have Pd(lb) as a NN display a substan- 
tial reduction in their chemical energy due to the large number of neighboring Pd atoms (the 
entire surface layer and the isolated Pd(lb) atom). Additional reductions in energy occur from the 
surrounding Ni atoms in the 2b and 3b planes. In spite of these contributions, however, there is a 
net total increase in the energy of the cell with respect to the case of a pure Ni (100) slab with a Pd 
surface layer (A E = 0.5530 eV). 
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That is not the case for the Pd(2b) atom shown in Fig. 3.b, since the number of Ni-Pd bonds is 
maximized. The Pd(S) atoms are barely affected by the Pd(2b) atom which in turn benefits greatly 
from the (maximized) number of Ni-Pd bonds thus created. The contribution of Pd(2b) is 0.2430 
eV, much lower than the Pd(lb) case (0.7697 eV) in spite of the fact that both atoms have almost 
identical coordination. In the Pd(2b) case, the chemical energy more than compensates for the 
near-bulk strain of the Pd(2b) atom. There is a net decrease in energy (-0.1205 eV) for this cell 
with respect to a pure Ni(100) slab with a Pd surface layer, thus favoring Pd atoms in this layer. It 
is also clear that, as a result of the competition between strain and chemical effects, an oscillatory 
profile emerges, as Pd atoms will populate alternating planes, in spite of the high strain, in order to 
maximize the gain in chemical energy. 

Figs. 3.c-d show the corresponding results for the (110) orientation. The most noticeable fea- 
ture is the negative contribution to AH of the Pd atom in the lb layer, as seen in Fig. 3.c. Added to 
the energy gains introduced by the affected surface Pd atoms, as well as the surrounding Ni atoms, 
a substantial lowering of the energy with respect to the case when the Pd(lb) atom is in the bulk is 
observed (-0.3087 eV). When the Pd atom is in the 2b layer (Fig. 3.d), the same behavior for the 
affected Ni atoms is observed (i.e., small negative contributions to the energy of formation). In 
this case, however, there is a substantial increase in the contribution of the Pd(2b) atom, thus lead- 
ing to a higher total energy for the cell than that obtained when this Pd atom is in the bulk (0.2087 
eV). As a result, Pd(lb/110) is favored, explaining the segregation of Pd to the top two layers for 
this orientation. 

3.2 Bulk alloy behavior 

The results from the previous section suggest that there is at least a weak tendency for order- 
ing in the Ni-Pd system. However, unlike many binary systems that exhibit a variety of ordered 
intermetallic phases, the Ni-Pd system is considered an ideal solid solution system [14]. However, 
experimental evidence for Ni-Pd does suggest the existence of at least short-range order effects 
[22-25]. Anomalous physical properties (electrical, magnetic, heat capacity, and other thermody- 
namic measurements) as a function of composition strongly hint at the existence of short-range 
ordering, particularly in compositions rich in Pd [23-25]. Moreover, a detailed study by Lin and 
Spruiell [22] suggests that local ordering occurs in Ni-50Pd and Ni-75Pd alloys after slow cooling 
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or aging at low temperatures (less than -673 K). These "ordered" regions consisted of clusters of 
like atoms in the {100} planes with roughly branched-rod shapes. Within the clusters, like atoms 
were linked primarily by second nearest neighbors. 

Since the BFS method is capable of performing bulk and surface analyses using the same 
parameter set, the modeling effort was extended to investigate the bulk structure of Ni-Pd alloys in 
order to further elucidate the possible role of short range order in controlling the near-surface 
structure of these alloys. For this purpose, large scale, temperature-dependent MC-NN simula- 
tions were performed for a variety of Ni-Pd compositions in order to shed light on the bulk struc- 
ture of such alloys. 

Fig. 4 contains interesting MC-NN results concerning the ordering behavior of several Ni-Pd 
alloys as a function of composition and temperature. Starting with a random initial state, repre- 
senting a high-temperature alloy, the temperature is lowered in 100 K steps, letting the system sta- 
bilize at any given temperature step. The most striking features of these results are the set of low 
temperatures obtained for the order-disorder transition for the different stoichiometries and the 
seemingly unstable ordering patterns observed. The low order-disorder transition temperatures are 
consistent with the general observation by Lin and Spruiell [22] that it is necessary to age alloys at 
very low temperatures, less than 675 K, before short range ordering was observed in Ni-50Pd and 
Ni-75Pd alloys. 

Further information on the degree and type of ordering can be obtained from the coordination 
matrices corresponding to the final states of the simulations (Table 2). The matrix element (by) 
of the coordination matrix a (b), for NN (NNN), respectively, denotes the probability that an atom 
of species i has an atom j as a NN (NNN). Consider, in particular, the Ni-75Pd alloy shown in Fig. 
4.b.The NN coordination matrix a resembles that for an LI 2 type ordered alloy. That is not the 
case for the NNN coordination matrix b, which shows a slight departure from the ideal values for 
LI 2 ordering. Visual inspection of the results indicates that this might not be a mere result of the 
particular simulation performed, but a consequence of the ordering of Pd atoms in rows along the 
[100] direction, with each Pd atom located at NNN distance from each other. In other words, the 
system does not appear to be evolving towards the tight ordering implied for an Ll 2 structure, but 
instead tends to form long strings of Pd atoms in specific directions. These simulations appear to 
be consistent with the x-ray study performed by Lin and Spruiell [22]. 

Fig. 4.c shows results of a simulation where the computational cell - representing a Ni-90Pd 
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alloy - was cooled down to T = 50 K. Not surprisingly, the high temperature results show mini- 
mum order, almost indistinguishable from a random distribution. As the temperature decreases, 
particularly below 600 K, Ni atoms start to group occupying NNN sites. It is only at room temper- 
ature, or below, that more noticeable changes occur, as the Ni atoms coalesce in what amounts to 
a small cluster of an Ll 2 -like NiPd 3 structure. 

The unusually low temperature for this transition can be traced to the nature of the Ni-Pd 
bonding and, from a modeling standpoint, can be understood by defining a small set of atomic 
configurations where a Ni atom occupies different possible sites relative to a fixed Ni chain 
embedded in a Pd alloy: (A) far away from the chain, (B-C) two possible NNN sites and (D-E) 
two possible NN sites, as defined in Fig. 5. a. After computing the energies of formation of these 
cells, an energy level spectrum (Fig. 5.b) can be built, highlighting the energy differences (relative 
to the lowest energy state) between the various configurations. The configurations with the lowest 
energy are those where the Ni atom joins the Ni chain in NNN sites, the seed of Ll 2 ordering, 
whereas the configurations with substantially higher energy are those where the Ni atom occupies 
NN sites, the seed for the formation of a Ni precipitate. The important feature is the proximity of 
the energy level of configuration A with respect to the other low-lying states, indicating that the 
bonding between Ni and Pd atoms is weak enough so as to barely differentiate states where Ni 
atoms follow a specific ordering pattern (configurations B and C) or are simply in solution (con- 
figuration A). Small thermal excitations can therefore break these bonds which are ultimately 
established at very low temperatures. 

Further understanding of this issue can be obtained from the atom-by-atom analysis of the two 
competing structures, as schematically shown in Fig. 6. In the first case (Fig. 6.a), two Ni atoms 
are separated by a distance large enough so as to ignore their mutual interaction. The local envi- 
ronment of each Ni atom is defined by the Ni atom itself and its 12 NN and 6 NNN. The contribu- 
tion of the Ni atom to AH is -0.0156 eV, and the contribution of each of the neighboring Pd atoms 
is -0.0020 eV/atom (with a negligible contribution from NNN Pd atoms). This results in a net con- 
tribution of -0.0803 eV for the two clusters. Both types of atoms have negative chemical energies. 
In the case of the neighboring Pd atoms, however, the chemical energy contribution per atom is 
small (-0.0020 eV). In the second case, shown in Fig. 6.b, where the two Ni atoms occupy NNN 
sites and the local environment now includes these two atoms and all their Pd NN and NNN. This 
situation is slightly less favorable for each of the two Ni atoms (-0.0150 eV) with respect to the 
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previous case in which they were separated by a large distance. The energy of the Pd NN atoms 
that are not shared by both Ni atoms remains unchanged. However, the four Pd atoms that are 
shared by both Ni atoms have now one additional Ni neighbor than the others, thus lowering the 
chemical energy from -0.0020 eV to -0.0085 eV, making the overall configuration more favorable 
than the separated Ni atoms in Fig. 6. a. The energy of formation of the cluster with two Ni atoms 
in NNN sites is -0.0967 eV, slightly lower than that for the two separate clusters (-0.0803 eV). 

There is one point in the previous example worth emphasizing. Although the Ni atoms display 
a large negative chemical energy contribution when surrounded by Pd atoms, this contribution is 
lessened as the number of Ni NNN grows. For example, a single Ni atom increases its energy (i.e., 
reduces its contribution to the energy of formation) by 0.0006 eV for each Ni NNN compared to 
an isolated Ni atom. While the energy contribution of Pd atoms favor increased coordination, 
inducing long range ordering, Ni atoms do not, instead preferring local ordering. If Ni atoms also 
benefited from this increased coordination, then the formation of an LI 2 structure would be more 
easily promoted, as both Ni and Pd atoms would tend to stabilize the structure, thus raising the 
ordering temperature. 

Taken as whole, the description of surface and bulk behavior by means of similar modeling 
tools helps one understand a variety of dissimilar issues in the framework of a simple, atom-by- 
atom approach. The behavior observed in the surface and near-surface region of Ni-Pd alloys can 
be rationalized in terms of these subtle interactions between Ni and Pd atoms. While the segrega- 
tion of Pd to the Ni surface can be basically explained by strain and surface energy concepts, other 
near-surface behaviors can only be explained by the additional fact that there is a slight attraction 
between Ni and Pd atoms. For example, the lb layer in a (110) surface is also heavily populated 
with Pd because the open characteristics of this surface orientation allow chemical effects to have 
a chance to lower the overall energy of the system. Likewise, the alternating Pd rich and deficient 
layers near the surface arise for the same reason. This type of behavior is usually associated with 
ordered alloys [15]. But, as shown in above, there is a weak tendency to order in these alloys since 
Pd atoms will lower their energy in the bulk by maximizing the number of Ni nearest-neighbors. 

4. SUMMARY AND CONCLUSIONS 
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A self-consistent model was presented for the description of surface and bulk properties of Ni-Pd 
alloys that not only reproduces existing experimental results but provides a coherent explanation for 
what was once a series of apparently conflicting data. It was shown that the competition between 
strain effects (leading to segregation of Pd to the surface) and chemical effects (leading to short-range 
ordering tendencies), must be taken into account to properly describe the top-layer-only Pd enrich- 
ment in (100) and (1 1 1) Ni-Pd surfaces vs. the two-layer Pd enrichment observed in a (1 10) termina- 
tion and the alternating enrichment and depletion of Pd in near-surface layers. Furthermore, the 
simulations and atomic analysis of bulk alloys show that the low short-range ordering temperature in 
bulk Ni-Pd alloys translates into a rapid smoothening of the oscillatory segregation profiles observed 
in the subsurface layers, especially with increasing temperature. 
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TABLE CAPTIONS 


Table 1: Experimental values for the lattice parameter, cohesive energy, and bulk modulus for the 
fee phases of Ni and Pd, and the resulting Equivalent Crystal Theory (ECT) parameters p, a, / and 
X. (see text and Ref. 17). 

Table 2: Nearest-neighbors and next-nearest-neighbors coordination matrices for Ni-50Pd, Ni-75 
Pd and Ni-90Pd alloys shown in Fig. 4 (T = 300 K). 

FIGURE CAPTIONS: 

Fig. 1: Segregation profiles for Ni-50Pd and Ni-5Pd alloys for the lowest index surfaces. Solid and 
dashed lines correspond to profiles obtained at 1000 K and 300 K, respectively. The horizontal 
axis denotes the plane (S = surface plane, lb = first plane below S, etc.) and the vertical axis indi- 
cates the concentration of Ni or Pd (in at.%). 

Fig. 2: a-c) Side view of a Ni (100) surface and d-f) a Ni (110) surface. Circles denote atoms in 
the plane of the page, while squares denote atoms one plane in and one plane out of the page. The 
Pd atom is indicated with a large black disk. Nearby Ni atoms (at NN or NNN distance from the 
Pd atom) are indicated with grey symbols. Other Ni atoms are indicated with open circles and 
squares. Each atom is labeled with its contribution to the total energy of formation of the slab, and 
the difference in energy with respect to a pure Ni(100) slab (in the corresponding tables). The 
total energy difference A E of the cell compared to a pure Ni surface is denoted for each cell. 

Fig. 3: Side view of a Ni (100) slab with a Pd surface layer and a Pd atom in the a) lb layer and b) 
2b layer, c) and d) show similar results for a Ni(l 10) surface. Circles denote atoms in the plane of 
the page, while squares denote atoms one plane in and one plane out of the page. Pd atoms are 
indicated with large black disks or squares. Nearby Ni atoms (at NN or NNN distance from Pd 
atoms) are indicated with grey symbols. Other Ni atoms are indicated with open circles and 
squares. Energy contributions of individual atoms to the energy of formation are indicated in the 
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corresponding tables. 


Fig. 4: Atomic distributions during a MC-NN simulation for Ni-50Pd, Ni-75Pd, and Ni-90Pd 
alloys. Starting from an initial, high temperature computational cell, the temperature is lowered to 
a final temperature of 300 K for Ni-50 Pd and Ni-75 Pd, and 50 K for Ni-90Pd. Equilibrium states 
at intermediate temperatures are shown, highlighting the range of temperature where the disorder- 
order transition takes place. Ni and Pd atoms are denoted with dark and open disks, respectively. 

Fig. 5: (a) (001) view of bulk configurations showing some possible locations of a Ni atom (grey 
solid squares) in the vicinity of a Ni-chain (black solid squares) in a Pd (open symbols) alloy: (A) 
away from the tip of the chain, (B-C) occupying a NNN site and (D-E) occupying a NN site, (b) 
Energy level spectrum showing the differences in energy of formation (in eV) between these con- 
figurations and the lowest energy state. 

Fig. 6: (a) Two Ni atoms (black squares) in a Pd (open symbols) alloy separated by a large dis- 
tance and their respective local environment (shaded area, containing the Pd NN to the central Ni 
atom), (b) Same, but with the two Ni atoms located in NNN sites. 
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Lattice 

Cohesive 

Bulk 






parameter 

(A) 

energy 

(eV) 

modulus 

(GPa) 

P 

a 

(A- 1 ) 

l 

(A) 

X 

(A) 

Ni 

3.524 

4.435 

148.81 

6 

3.015 


0.7587 

Pd 

3.89 

3.94 

156.23 

8 

3.612 


0.666 


Table 1 
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Nearest- 

neighbors 

Next- 

nearest- 

neighbors 



Ni 

Pd 

Ni 

Pd 

Ni-50Pd 

Ni 

0.33 

0.67 

0.83 
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Pd 

0.67 

0.33 

0.17 
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Ni-75Pd 

Ni 
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0.73 

0.27 


Pd 

0.33 

0.67 

0.09 

0.91 

Ni-90Pd 

Ni 

0 
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0.20 

0.80 


Pd 

0.11 

0.89 

0.09 

0.91 


Table 2 
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Ni-50Pd (100) 



surface (S) 


Ni-5Pd (100) 



Number of layers below the 
surface (S) 


Ni-50Pd (110) 



Number of layers below the 
surface (S) 


Ni-5Pd (110) 



Number of layers below the 
surface (S) 


Ni-50Pd (111) 



Number of layers below the 
surface (S) 


Ni-5Pd (111) 



Number of layers below the 
surface (S) 


Fig. 1 
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